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EFFECT OF BLADE -THICKNESS TAPER ON AXIAL -VELOCITY DISTRIBUTION 
AT THE LEADING EDGE OF AN ENTRANCE ROTOR -BLADE ROW WITH AXIAL 
INLET, AND THE INFLUENCE OF THIS DISTRIBUTION ON ALINEMENT 
OF THE ROTOR BLADE FOR ZERO ANGLE OF ATTACK 
By John D. Stanitz 

SUMMARY 

A method Is developed for estimating the effect of blade -thickness 
taper on the inlet axial -velocity distribution of an entrance rotor- 
blade row with axial inlet, and the influence of this velocity dis- 
tribution on the alinement of the rotor blade for zero effective angle 
of attack (that is, zero blade loading at the nose). This alinement 
of the blade requires a deviation between the angle of the blade cam- 
ber line direction at the inlet and the upstream relative flow direc- 
tion. The method is developed for compressible and incompressible non- 
viseous fluids, and results are presented for incompressible flow into 
a plane, two-dimensional cascade and for compressible flow into an 
entrance rotor -blade row with tapered blades. It is concluded that, for 
the entrance rotor -blade row investigated, blade taper has a large effect 
on the inlet deviation angle; whereas compressibility has a small effect, 
except perhaps at the hub, and the Tip stream relative flow direction also 
has a small effect. 


INTRODUCTION 

If the blade-element profile is to be set at the desired angle of 
attack at each radius in any blade row of a turbomachine, it is first 
necessary to know the blade-element alinement for zero "effective" angle 
of attack, that is, for zero loading at the inlet. In many cases this 
alinement is considerably different from the upstream relative flow 
direction, and if It is not known, improper angles of attack may result. 
Thus, rapid blade loading, accompanied by high peak velocities near the 
nose of the blade-element profile, may occur. The subsequent decelera- 
tion from these high velocities on the blade surface can result in local 
boundary-layer separation and accompanying mixing losses. For a blade 
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raw with, an essentially hlimt leading edge, the gas passes from the 
"approach" station, which is immediately upstream of the blade leading 
edge, to the "inlet" station, which is immediately downstream of the 
blade leading edge and has finite "effective" blade thickness. Aline- 
rnent of the blade-element profile for zero angle of attack, that is, 
for essentially zero loading at the inlet station, trill depend on the 
inlet axial velocity (at the inlet station) and must be such that the 
whirl (product of absolute tangential velocity and radius) is unchanged 
from its upstream value along the same streamline. For "entrance" rotor- 
blade rows, which are not preceded by guide vanes and which are the sub- 
ject of this report, the upstream whirl is zero. For this type of rotor- 
blade row, which occurs in centrifugal- and axial-flow compressors that 
do not have inlet guide vanes, the zero angle of attack condition is 
achieved when the blade profile at each radius is alined to preserve 
zero whirl at the inlet station immediately downstream of the leading 
edge. 


In the design of entrance rotor -blade rows with axial inlets it is 
common practice, although usually incorrect, to assume a uniform approach 
axidl -velocity distribution. Actually there are a number of factors that 
can result in nonuniform distributions. Among these factors are: 

(l) boundary layer on the hub and casing, (2) meridional curvature of the 
hub and casing -upstream and/or downstream of the blade leading edge, 

(3) blade curvature, (4) radial variations in angle of attack, and 
(5) blade taper (radial variation in blade thickness). 

The influence of boundary layer on the approach axial -velocity dis- 
tribution is obvious; likewise, meridional curvatures of the hub and 
casing should be expected to cause nonuniform distributions in which 
the highest velocities occur at the walls with larger convex curvature. 
The effect of hub and casing curvatures downstream of the blade leading 
edge, such as occur in centrifugal and mixed-flow compressors, has been 
noted in reference 1. Blade curvature and radial variations in angle of 
attack can result in radial displacements of the meridional stre aml ines, 
which induce nonuniform approach axial-velocity distributions. Likewise, 
blade taper requires radial displacements of the flow and therefore 
induces variations in the approach axial-velocity distribution. 

It will be assumed that the approach axial -velocity variations, and 
therefore the inlet axial-velocity variations, induced by these various 
factors can be treated separately, and that the resulting distribution 
can be determined with reasonable accuracy by a simple superposition of 
the several variations. It is, then, the purpose of this report to 
investigate, a n alytic al l y , the variation in inlet axial-velocity distri- 
bution due to blade taper alone, and to indicate the influence of this 
distribution on the alinement of the rotor blade -element profile for 
zero angle of attack at each radius. The investigation includes both 
compressible and Incompressible fluids; for compressible fluids the 
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condition of sonic relative velocity at the inlet station is also con- 
sidered. The object of this investigation, made at the NACA Lewis lab- 
oratory, is to obtain approximate values of the phenomena involved and 
improved knowledge of the physical mechanisms. 


METHOD OF ANALYSIS 

A method is developed for estimating the variation in inlet axial- 
velocity distribution induced by the blade taper in an entrance rotor- 
blade row. In addition, a method is developed for estimating correc- 
tions to the rotor -blade inlet angles so that the blade is alined with 
the flow, that is, so that the blade loading is approximately zero at 
the Inlet. 


Preliminary Considerations 

Statement of problem . - If all factors except blade taper affecting 
the inlet axial -velocity distribution in an entrance rotor -blade row are 
eliminated, the problem of the inlet axial -velocity distribution reduces 
to axial flow through an annulus of concentric circular cylinders in 
which Is located a cascade of straight blades with specified taper 
lying, In effect (see appendix B), in the meridional (axial -radial) 
plane and extending downstream to infinity. Such an annulus is shown in 
figure 1 together with a meridional streamline that is displaced radially 
outward by the blade taper. (The symbols used in fig. 1 and elsewhere in 
this report are defined in appendix A. ) Par upstream and downstream of 
the blade leading edge the streamlines are parallel to the axis and the 
axial velocity is uniform from hub to casing. In the vicinity of the 
leading edge, however, the radial displacement of the streamlines induces 
a variation in axial -velocity distribution from hub to casing at both the 
approach and the inlet stations. This variation in velocity distribution 
depends on, and will be estimated for, specified values of the uniform 
upstream velocity, blade taper, and hub-tip ratio. The variation in 
inlet axial -velocity distribution requires an adjustment from hub to tip 
in the rotor -blade inlet a n gle (inlet station, fig. l) in order to 

achieve zero blade loading at the inlet. 

Assumptions . - As already stated, it is assumed that the effects of 
blad e taper on the inlet axial -velocity distribution can be treated 
(with reasonable accuracy) separately from the effects of other factors, 
and that the resultant effect of ad 1 factors can be determined by simple 
superposition. (However, only the effect of blade taper is considered 
herein.) In addition, the flow is assumed to be axially symmetric. 

Also, the leading edge of the blade is assumed essentially blunt 
(fig. l) , although shaped to avoid entrance losses. Both compressible 
(subsonic) and inconpressible nonvisous fluids are considered. Further 
assunptions are introduced as needed. 



4 


NACA TN 2986 


Outline of method . - The variation in inlet axial -velocity distri- 
bution that results from blade taper in an entrance rot or -blade row is 
obtained from a stream function distribution at the inlet station that 
is related,, by assumption, to the known upstream and downstream distri- 
butions. (The validity of this assumption is checked in appendix B.) 
After the inlet axial -velocity distribution from hub to casing is known, 
the inlet deviation angle (inlet station, fig. l) of the rotor -blade 

inlet angle p ± from the upstream relative flow direction P u is deter- 
mined in such a way that the absolute tangential, velocity is unchanged 
(zero), that is, so that the blade is not loaded at the inlet. It is 
supposed that knowledge of this deviation angle is needed, as a refer- 
ence, in order to set the blade -element profile at the desired angle of 
attack. 

CoorrHrtntes and velocity components . - A point between the hub and 
casing in the meridional plane is designated by its r,z coordinates, 
where r is measured from the axis of the annulus and z is measured 
from the leading edge of the blade, positive in the direction of flow 
(fig. l) . These coordinates are dimensionless, being expressed as 
ratios of the casing radius (so that r is 1.0 at the blade tip, for 
exanple) . The hub radius is designated by r^. 

The velocity components in the r and z directions are v and 
u, respectively. These velocity components are dimensionless, being 
expressed for incompressible flow as ratios of the upstream velocity 
(so that for incompressible flow is equal to 1.0, for example), 

and being expressed for compressible flow as ratios of the upstream 
stagnation speed of sound. 

The angular velocity co of the rotor is dimensionless, being 
expressed in such a way that the dimensionless wheel speed cur is a 
ratio of the upstream velocity for incompressible flow and a ratio of 
the upstream stagnation speed of sound for compressible flow. 

As a result of the wheel speed, the flow relative to the rotor has 
the direction p measured from the axial direction on a developed 
cylindrical surface (fig. l) . The difference between the upstream rela- 
tive flow direction (3 U and the relative flow direction p is defined 
as the deviation angle e (fig. l) . 

€ = Pu - P (1) 

Blade taper . - The blade taper is determined 'by the variation in 
the passage width ratio S, which is defined at each radius as the ratio 
of the circumferential passage width between adjacent blades to the 
blade spacing (fig. l) . Thus, if a is the angular spacing of the 
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blades about the z-axis, 5rn is the passage width, and the circum- 
ferential blade thickness t is given by (fig. l) 


t = (l-S) ra 


( 2 ) 


Thus, if the blade taper is defined as - (dt/dr), equation (2) results 
in 


dt 

dr 





(3) 


In this report the passage width ratio 5 will be specified as a func- 
tion of r by the power series 

5 = a + br + cr 2 + . . . (4) 


Inlet Axial -Velocity Distribution 


Stream function ¥ . - As a result of continuity considerations a 
stream function W is defined (appendix B) such that 


d¥ _ Srpu 
<5r ~~ A >Jr 


(5) 


where for isentropic flow the density p is given by 

1 

P - (l - Eji (6) 


and where, in terms of upstream flow conditions, 

■** - <vv (-r~) 

or, in terms of downstream flow conditions. 


(7a) 


= Pd% |[d- r h 2 ) + I d-^h 3 ) + f d- r h 4 ) + • • •] (7t) 

The stream function Uf has been defined (by the introduction of Aijr) so 
that it varies from 0 to 1.0 between the hub and the casing, respectively. 
The density p Is dimensionless, being expressed as a ratio of the 
upstream stagnation density. For incompressible flow p is 1.0. 
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Upstream distribution of iff . - Far upstream of the blade inlet, 
both p and u are constant and 6 is 1.0; so that equation (5) can 
be integrated to give, when combined with equation (7a), 




(8) 


Equation (8) gives the upstream distribution of If as a function of r. 

Downstream distribution of W . - Far downstream of the blade inlet, 
both — p and u are constant (because all factors other than blade 
taper have been neglected) and 8 is given by equation (4); so that 
equation (5) can be integrated to give, when combined with equation (7b), 



(9) 


Equation (9) gives the downstream distribution of ¥ as a function of r. 

Inlet distribution of ¥ . - At the blade inlet, both p and u are 
unknown and variable and 6 changes suddenly from an upstream value of 
1.0 to its prescribed downstream values . Under these unknown conditions 
equation (5) cannot be integrated; instead it is assumed that the inlet 
value of ¥ at each radius is equal to the average of its upstream and 
downstream values at the same radius. 

- \ % + V (10) 

This basic assumption is shown to be valid by a series of relaxation 
solutions in appendix B. Because the density is constant at both the 
upstream and downstream stations, the distributions of ¥ u and 

are independent of compressibility effects (see eqs. (8) and (9)). 

Equation (10) therefore assumes that is also independent of com- 

pressibility effects. Inasmuch as the density is not constant at the 
inlet station, this assumption is also investigated in appendix B. 

Inlet axial-velocity distribution . - The inlet axial -velocity dis- 
tribution is obtained from the radial derivative of the distribution of 
¥j_ given by equation (10). Thus, from this derivative and equa- 
tions (5), (7a), (8), and (9), 
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Pi u i = 


Pu^u 


§ + - 
o a 




(l-r h 2 ) + | (l-r h 3 ) + | (l-r h 4 ) 


(H) 


The distribution of u^ from hub to casing is given by equation (ll) in 
which is related to u^ by equation (6) and 5 is a specified 

function of r given by equation (4). 


Rot or -Blade Alinement 

A method is developed for computing the radial variation in inlet 
deviation angle e i required to avoid loading of entrance rotor blades 

at the inlet. Although the method, as developed, takes into account the 
effects of blade taper only, it can be extended, by methods to be dis- 
cussed later, to include the effects of other factors. . 

Inlet deviation angle e-j_. - For an entrance rotor blade the rela- 
tive upstream flow direction p u is given by (fig. l) 

tan p u = - 1 (12) 

"u 

(in order to avoid negative values of |3, osr is considered positive as 
shown in fig. 1.) Equation (12) can also be expressed as 

tan P u = r(tan P u ) 1>0 (l2a) 

where the subscript 1.0 refers to the casing at which r is 1.0. At the 
blade inlet, immediately downstream of the leading edge, 

i q GOT 

tan pj = — 

1 u i 


which, frorfL equation (12), becomes 


tan p i = — tan p u 


Also, from equation (l), 


tan p u - tan 

l 811 i + tan P u tan 0^ 


( 13 ) 
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which, from equation (13), becomes 


tan e i 


(l . ^ ta n p u 

^U. o 

1 + — tan^ p„ 

Ui u 


(14) 


The distribution of inlet deviation angle €± from hub to casing is 
given by equation (.14), in which u^ is given by equation (ll) and 0 U 
is given by equation (12a). 

Special case. - If the flow is plane and incompressible, r h and p 
are 1.0, and equation (ll) reduces to 



so that equation (14) reduces to 


(1 - 5) tan 

tan = ■= (15) 

1 + 6 tan 

Equation (l5) has been developed by Weinig (ref. 2, p. 120) . 

Condition for sonic inlet relative velocity . - It can be shown that, 
if the blades are alined for zero loading at the nose, the resultant 
inlet relative velocity is sonic when 

t \ 2 y +1 2 _ 

(cur) + -- uj_ = 1 

which, from equation (12), becomes 

u t 2 = (1 - u u 2 tan 2 p u ) (16) 

From equations (6), (ll), and (16), 
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IT+ 1 ^2 


- tan 2 p u 


1 ' m (1 " ^ tan2 Pu) 


1 - £— - -n 2 

J- 2 ^ 


r-i 


2 

5 + 


1 - r> 


| (l-r^ 2 ) + ~ (l-r h 3 ) + j (l-r h 4 ) + . . . 


(17) 


oj 

i 


6 


Equation (17) is the condition for which the resultant relative inlet 
velocity is sonic. This condition is of interest, as a limit, if the 
designer wishes to avoid supersonic inlet relative velocities. 


Numerical Procedure 

Specified conditions . - In order to investigate the effect of "blade 
taper on the inlet axial -velocity distribution of an entrance rotor- 
blade row and the influence of this distribution on the inlet deviation 
angle c^, the following conditions must be specified: 

< tan Ml.O I 


r h 

5 = 5(r) j 


(18) 


For compressible flow is expressed as a ratio of the upstream stag- 

nation speed of sound; for incompressible flow ^ is 1.0. The speci- 
fied distribution of 8 as a function of r is given by the coeffi- 
cients a, b, c, . . . in equation (4). The specified distribution of 
8 is related to the blade taper - (dt/dr) by equation (3) . 

Calculations . - The following procedure is used in the calculations 

(l) If the? fluid is compressible, and if the designer wishes to 
avoid supersonic inlet relative velocities, the conditions (u^ and p u ) 
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for which the inlet relative velocity is sonic should he determined 
from equation (17) . 

(2) If from step (l) the assigned value of u u does not result in 
sonic relative velocity, compute u^ from equation (ll) in which p is 
related to u "by equation (6). If the fluid is incompressible, p is 

1 . 0 . 

(3) Compute from equation (14) . 

If other factors besides blade taper affect u^ and e^, it is 
first necessary to compute the approach axial-velocity distribution Uq, 
(at the approach station, fig. l) resulting from each of these factors 
considered separately. For example, the approach axial-velocity distri- 
bution Ug resulting from blade taper alone is given by the continuity 
equation 

Pa^a = 5 Pi u i ( 19 ) 

where p-jU^ is given by equation (ll) and p a is related to Ug by 
equation (6). The variations in Ug from that result from each of 

the factors separately are then added to obtain the resultant distribu- 
tion of Ug. The final distribution of u^ is then obtained from equa- 
tion (19), using the resultant distribution of Ug, and the distribution 
of inlet deviation angle is then determined by equation (14). 


RESULTS AND DISCUSSION 

Results are presented for incompressible flow into a plane, two- 
dimensional cascade and for compressible flow into an entrance rotor - 
blade row with tapered blades. 

Plane cascade with incompressible flow . - For incompressible flow 
into a plane, two -dimensional cascade, equation (15), developed by 
Weinig, relates the inlet deviation angle € ± to the upstream flow 
direction flu and the passage width ratio 5. The relation is plotted 
in figure 2. As expected, c i increases with decreasing values -of 8. 

It is interesting to note that for values of 8 commonly occurring in 
practice, say 0.7 or greater, the values of are 10° or less. For 

values of 5 as low as 0.5, such as might occur near a blade root, the 
values of can be as high as 20°. 
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Entrance rotor-blade row -with Hade taper . - The more general case 
of compressible flow into an entrance rotor -"blade row with the same hub- 
tip ratio and blade taper u§ed in example V of appendix B is considered. 
First, in order to avoid specifying upstream flow conditions (uu and 
(p u )l q) that result in local sonic flow at the inlet station, the criti- 
cal value of Uu (at each radius) for which the relative inlet velocity 
is sonic, determined by equation (17), is plotted in figure 3 as a func- 
tion of r for three values of (Pu)i.O* r ^ ie variation in this critical 

value of u u with r in figure 3 is not large for this thickness dis- 
tribution, because the blade speed increases the inlet relative velocity 
toward the casing and the blade taper increases the velocity toward the 
hub. (This taper effect must be large because, as will be shown, 

increases toward the hub and this increase tends to reduce the inlet 
relative velocity there.) For (Pu)i.o e l ua l do 60° the variation in 

u u with r lies between 0.45 and 0-50, indicating that sonic relative 
velocity occurs almost simultaneously from hub to casing. For (P u )]_,o 

equal to 50°, sonic inlet relative velocity first appears at the hub 
(as indicated by the minimum value of u u at this radius) } whereas 
for (p ), n equal to 70°, sonic velocity first appears at the casing. 

The variation in inlet deviation angle with radius for various 

values of u u and (Pu)i.o was determined by equations (ll) and (14) 
and is plotted in figure 4. In figure 4(a) the effect of on 

is shown for (P u )qQ S'! 1181 ! to 60°. A maximum value of 0.4 for 

was chosen in order to avoid the sonic inlet relative velocity pre- 
dicted by figure 3. The higher values of toward the hub indicate 

a large effect of blade taper on The spread in the curves for var- 

ious values of u^ indicates that, for the example investigated, the 
effect of compressibility on £j_ is not important, except perhaps near 

the hub (where, however, the relaxation solutions in appendix B show 
that the estimated values of Uq (and therefore e q) are too high, 
especially for high subsonic values of u^). 

In figure 4(b) the effect of (P-a^.o on € i is s ^ own for ^ 

equal to 0 (corresponding to incompressible flow). Again the curves 
indicate a large effect of blade taper on eq* The spread in the curves 

for various values of (P u )q q dedicates that for the example investi- 
gated the effect of (P u )q q on e i is n °d important. This result 
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agrees with that for plane cascades ■with incompressible flow where 
(fig. 2), for the range of 8 considered in figure 4(h), only small 
variations in occur for large variations in f3 u . However, as the 

value of 8 decreases it is observed in figure 2 that the influence of 
3 on £i becomes more pronounced. 

Further considerations . - In the practical application of the 
methods developed in this report for estimating the radial variation in 
€j_, certain difficulties may arise with regard to the selection of 8 

as a function of r. These diff iculties arise if the blade nose is 
faired (fig. 5) , instead of being blunt (fig. l) , so that 8 is a 
function of z as well as of r. In this case it is necessary to 
assume some "effective" variation in 8 with r, which might, for 
example, be the variation 2 or 3 percent of the blade height downstream 
of the blade leading edge. Another difficulty, discussed in appendix B, 
is the breakdown for high values of u^ in the basic assumption that u^ 

can be computed from the distribution of given by equation (10). 

In both these cases, however, the estimated values of u^, and 

can serve as a guide to the designer (if properly interpreted). 

In this analysis it has been assumed that the blade -element pro- 
files have blunt noses (or, in the event that the noses are not blunt, 
that some "effective" blade thickness can be assumed) . It is therefore 
of interest to examine the inlet profiles that Welnig (ref. 2) has 
developed for the type of two-dimensional cascade blades 'considered in 
figure 2. These profiles have been developed to give a continually 
accelerating flow from the stagnation point at the nose to the maxi Timm 
velocity in the downstream channel between blades. Such a velocity 
distribution eliminates the possibility of local bound ary- layer separa- 
tion near the nose. From equations on page 122 of reference 2, the pro- 
file coordinates, x along the blade chord and y normal to the chord, 
are given by 

X = 8 C °« Se ° [ (l-S)los Pl - P u] Pl ' fc] tan <v} 

( 20 ) 

where, for the type of cascade being considered, 

tan = 8 tan P u (2l) 

and where the x and y coordinates are expressed as ratios of the 
circumferential blade spacing. Examples of the blade profiles are given 
in figure 5 for three Values of 8 and for P u equal to 45°. The 
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radius of curvature R at the "blade nose, expressed as a ratio of the 
circumferential "blade spacing, is given "by (ref. 2, p. 12l) 

8(l-8)^ cos P-r 

R = s—— ( 2 ; 

^(5^ + tan^ p^) 

This radius might he of some help in approximating the desired blade - 
element profile, and it is therefore plotted in figure 6 as a function 
of .8 for various values of P u . 


SUMMARY OF RESULTS AND CONCLUSIONS 

A method is developed for estimating the effect of "blade taper on 
the inlet axial -velocity distribution of an entrance rotor -blade row 
with axial inlet and the influence of this velocity distribution on the 
alinement of the rotor blade for zero loading at the nose. This aline- 
ment of the blade results in an inlet deviation angle of the blade- 

camber-line direction at the inlet from the upstream relative flow direc 
tion. The method applies to both compressible and incompressible nonvis 
cous flows, and an expression is developed for predicting local sonic 
relative velocity at the inlet. The method is applied to plane, two- 
dimensional cascades with incompressible flow and to the more general 
problem of compressible flow into an entrance rotor-blade row with blade 
taper. It is concluded that for the entrance rotor-blade row investi- 
gated: 

1. Blade taper has a large effect on 

2. The effect of compressibility on is small, except near the 

hub . 

3. The effect of upstream relative flow direction on is small. 

4. The upstream velocity for which the relative inlet velocity is 
sonic at a given radius does not vary greatly from hub to casing, (in 
one example sonic relative velocity occurred almost simultaneously from 
hub to casing.) 


Lewis Flight Propulsion Laboratory 

National. Advisory Co mmi ttee for Aeronautics 
Cleveland, - Ohio, May "19, 1953 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: (All symbols are 

dimensionless; velocities are expressed as ratios of the -upstream abso- 
lute stagnation speed of sound for compressible flow , and as ratios of 
the upstream axial velocity for inconpressible flow; distances are 
expressed as ratios of the casing radius , unless otherwise specified. ) 

a,b,c, ... coeff icients, prescribed to define 6 by eq. (4) 

R radius at nose of two- dim ensional blade profile, expressed 

as ratio of blade spacing in two-dimensional cascade , 
eq. (22) 

r radius, measured from axis of annulus, fig. 1 

t circumferential blade thickness, eq. (2) 

u,v axial and radial velocities, respectively 

x,y coordinates of two -dimensional blade profile, measured 

along blade chord and normal to it, respectively, and 
expressed as a ratio of blade spacing in two-dimensional 
cascade, eq. (20) 

z axial distance, measured from leading edge of blade, posi- 

tive in direction of flow, fig. 1 

P relative flow angle, fig. 1 

X ratio of specific heats 

6 passage width ratio, ratio of circumferential, passage width 

to circumferential blade spacing at same radius 

£ deviation angle, eq. (l) 

p density, expressed as ratio of upstream stagnation density, 

eq. (6) 

a angular spacing of blades around z-axis, fig. 1 

^ stream function, eq. (B2) 

/hjr constant, defined by eq. (B3) 
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co angular velocity of rotor , expressed in such a way that the 

dimensionless wheel speed cor is a ratio of upstream 
velocity for incompressible flow, and a ratio of upstream 
stagnation speed of sound for compressible flow 


Subscripts: 

a approach, immediately upstream of blade leading edge (fig. l) 

d .downstream, far downstream of blade leading edge 

h hub 

i inlet. Imme diately downstream of blade leading edge (fig. l) 

u upstream, far upstream of blade leading edge 

casing, where r is 1.0 


1.0 
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APPENDIX B 

RELAXATION SOLUTIONS FOR DISTRIBUTION OF ¥ 

The validity of the "basic assumption that, for "both compressible 
and incompressible flows. 


= I * %) 


( 10 ) 


is investigated "by comparing the results of this assumption with the 
results of seven relaxation solutions for flow thro ugh annuli of the 
type shown in figure 1. 


Method of Analysis 

If it is assumed that the rotor "blade profiles at the inlet, and, 
for purposes of this investigation, everywhere downstream of the "blade 
leading edge, can he alined with the flow in such a manner that no work 
is done on the fluid, then the fluid motion in the meridional (axial- 
radial) plane is irrotational, so that 


chi chr „ 

SF ~ Si" 0 


(Bl) 


Also, for axially symmetric flow, which assumes that the number of Blades 
becomes infinite but that the blockage effect of the blades (determined 
by 8) is unchanged, the continuity equation becomes 

(5rpv) + ^ (Srpu) = 0 


This continuity equation is satisfied by a stream function ¥, which is 
defined by 


_ Srpu 
cEr ~ Zty 

_ Srpv 
cEz ~ ~ Zbjr > 


(B2) 


where the density p is related to u and v by the usual isentropic 
equation 
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= £l - (u2 + v2) 


and where the constant Ai{r is defined by 

v 1.0 


A = 


6rpu dr 


(B3) 


so that the stream function UT varies from 0 at the hub to 1.0 at the 
casing. In terms of upstream flow conditions, where p and u are 
constant and 5 is 1.0, equation (B3) can be integrated to give 


At = P u u^ 



(7a) 


Also, in terms of downstream flow conditions, where p and u are con- 
stant and 5 is specified by equation (4), equation (B3) can be inte- 
grated to give 


^ ~ Pd 1 ^ \_2 ( 1-r h 2 ) + 3 ( 1 “ r h 3 ^ + 4 ( 1-r h 4 ^ + • • -J 
Equations (Bl) and (B2) combine to give 


(7b) 


+ W d In (5rp) _ 3¥ a In p 

^.2 ^ z 2 Sr ot Sz 


(B4) 


Equation (B4) determines the distribution of the stream function ¥ from 
hub to shroud, and the resulting velocity distribution is determined by 
equations (B2). It is noted that the form of equation (B4) is indepen- 
dent of the flow direction (or blade angle) (3 so that, in effect, the 
rotor blades can be considered to lie in the meridional plane, in which 
case they must also be stationary. That is, if the rotor blades are 
ali n ed to do no work on the fluid, they have no effect on the absolute 
meridional flow other than that due to blade taper, which effect is the 
same for stationary blades lying in the meridional plane. 


Numerical Procedure 

Equation (B4) was solved by relaxation methods (ref. 3) for speci- 
fied geometric and flow conditions. A square grid with 0.0625 spacing 
was used, and the finite-difference equations were based on a second- 
degree polynomial for ¥. The residuals were reduced to give five-figure 
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accuracy for the value of ¥ at each grid point. The boundary values 
of ¥ were 0 and 1.0 along the hub and casing, respectively, and the 
relaxation solutions were extended upstream and downstream of the lead- 
ing edge until flow conditions became uniform. 


Numerical Results 

Seven relaxation solutions were obtained for the conditions speci- 
fied in the following table: io 


Example 

^u 

B 

6 

Comment* 

I 

0 

( incompressible ) 

0.625 

-0.25 + 1.2r 

\ = 0- 5 

II 

0 

( incompre s sible ) 

.625 

0.0167 + 0. 9333r 

\ = 0.6 

III 

0 

( incompressible ) 

.625 

0.2833 + 0. 6667r 

Sh = 0.7 

IV 

0 

(incompressible) 

.500 

0.05 + 0.9r 

\ = 0.5 

V 

0 

(incompressible) 

.625 

-1.4556+4. 7611r-2.3556r 2 

8^0. 6 =0 
^ \ dr /1.0 

VI 

0.44 

.625 

-1 . 4556+4 ;7 611r -2 . 3556r 2 

^°- 6 '(i)i.o =0 

VII 

0.52 

.625 

-1 . 4556+4 . 761hr -r2 . 3556r 2 

6 h= 0 ' 6 '(i) 1 .o =0 
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u u (0, 0.44, and 0.52). The discontinuity in velocity distribution at 
the leading edge results from the sudden change in 5. This discon- 
tinuity becomes noticeably more severe as the value of u^ increases. 

Also, the large axial gradients of u downstream of the leading edge 
are an indication of large axial gradients of e, which must exist if 
the blade is to remain unloaded. The largest value of u^, 0.52, is 

the maximum value for which the relaxation solution would converge and 
corresponds to the condition of sonic velocity at the blade root. 
(Actually, eq. (17) with equal to zero indicates a value of 0.4500 

for sonic velocity on the hub at the blade leading edge. However, 
because grid points were not located on the blade leading edge, sonic 
velocity was not encountered in the relaxation solutions until u u was 
greater than 0.52.) 

Flow conditions along the blade leading edge as obtained from the 
relaxation solutions (I through VII ) are compared In figures 9 and 10 
with the estimated flow conditions based on the assumed relation given 
by equation (10). In figure 9, the estimated distribution of ¥q 

given by equation (10) is compared with the distribution of '<?£ obtained 

by the relaxation solutions. The agreement is seen to be excellent. 
However, it is the accuracy of the radial derivative of ¥j_ , from which 
Uj[ is obtained, that is of interest. The upstream value of this deriva- 
tive at a given radius divided by the inlet value at the same radius is 
equal to p^u^Bp-jU^. The estimated distribution of this quantity, 

obtained from equation (ll), is compared in figure 10 with the distribu- 
tion obtained by the relaxation solutions. Again the agreement is sat- 
isfactory, except perhaps for large values of p^u^Sp^u^, which occur 

near the hub. For these values of p u u u /5p^u^ the estimated value of 
PjU^ is too large so that u^ is too high, especially for high sub- 
sonic values of u^ where, because PqUq is almost constant, very 
large cha n ges in u^ are required to achieve small changes in PqUq* 
However, except for high subsonic values of uq near the hub, it is 

concluded that the basic assumption represented by equation (10) is 
valid for both compressible (subsonic) and incompressible flow. 
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Figure 3. - Variation in u u required at radius r for sonic velocity 
in inlet. Example Y, equation (17). 
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(a) (Pu)l.O> 60 °* 

Figure 4. - Variation in inlet deviation angle e i with radius r for 
example V with various values of (P u )j_ 0 and u^. Equation (14). 
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to 

r- 

05 

W 


I 



.625 .750 .875 1.000 

r 

(b) 0 (corresponding to incompressible f lov) . 

Figure 4. - Concluded. Variation in inlet deviation angle with radius 

r for example V with various values of (P u ) 1 q and u u * Equation (14). 




Figure 6. - Variation in nose radius R with passage ratio 
5 for various relative upstream flow angles £ u . Plane, 

two-dimensional, incompressible flow) equation (22). 
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3tr*am funotloo, Blad* 

7 lcadlcf 



Figure 7. - Concluded. Relaxation solutions for stream line distribution upstream and 
downstream ctf leading edge of entrance rotor-blade row with blade taper. Dashed lines 
show radial spacing of streamlines far upstream and far downstream of leading edge. 
5 1.0> °- 95 ; incompressible flow. 








(o) ra. Uu ■ 

». - % » •«-““• “» 
. r,«. n.. 0.80. ft.. 0.625 1 6 - -1.8656 +■ 4.7611P - E.SWWr 2 , o«pM».iM« flow. 
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Figure 9. - Comparison of stream function obtained along 

leading edge of blade by relaxation solution with approxi- 
mate stream function Wi given by equation (10) . Values 
of r ± : 0.625* 0.6875, 0.750, 0.8125, 0.875, 0.9375, and 
1.000$ in addition, 0.500 and 0.5625 for example IV. 
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Example 



.8 .9 1.0 1.1 1.2 1.3 

Relaxation value of Pu u u/( & Pi u i) 

Figure 10. - Comparison of p u u u /(Sp i u i ) obtained along leading 
edge of blade by relaxation solution with Pu u u/^Pi u i) given 
by equation (ll). Values of r^: 0.625, 0,6875, 0.750, 0.8125, 
0.875, 0.9375, and 1.000$ in addition, 0.500 and 0.5625 for 
example IV. 
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